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ABSTRACT 
We  investigate  the  corrections  to  the  representation 
of  the  joint  distribution  of  q+j£     particles,  nQ+rf>     b7 
the  product  n  n»      for  large  separation  between  the  sets 
of  q  and  s&     particles.   For  a  system  in  which  there 
exists  a  "finite  correlation  length,"  we  find  explicitly 
the   l/N   correction  term  to  the  simple  product,  where   N 
is  the  number  of  particles  in  our  system.   When  q+A     is 
equal  to  two,  this  expression  reduces  to  that  familiar  from 
the  Ornstein-Zernike  relations  for  scattering  of  light  from 
a  fluid.   In  a  uniform  gas,  our  derivation  also  yields  the 
complete   l/N  dependence  of  equilibrium  distributions. 
Our  result  on  the  asymptotic  form  is  then  used  to  determine 
the  low  order  distribution  functions  for  an  equilibrium 
system  of  varying  density,  as  well  as  for  a  non-equilibrium 
system  represented  by  a  local -equilibrium  ensemble.   These 
distribution  functions  are  shown  to  be  governed  by  the 
temperature  and  density  in  the  vicinity  of  the  molecules 
considered.   We  find  as  expected  that  the  two-body  distri- 
bution function  coincides,  to  within  quadratic  terms  in 
the  gradients,  with  its  equilibrium  value  for  a  uniform 
system  at  the  temperature  and  density  of  the  midpoint. 
For  the  higher  order  distributions,  correction  terms  linear 
in  the  gradients  are  found. 
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LONG  RANGE  CORRELATIONS  IN  A  CLOSED  SYSTEM  WITH 

APPLICATIONS  TO  NONUNIFORM  FLUIDS 

lo  Introduction. 

One  of  the  fundamental  concepts  of  macroscopic  physics  is 

that  of  a  homogeneous  system.   The  state  of  such  a  system  is 

completely  described  by  a  set  of  intensive  parameters  which  make 

no  reference  at  all  to  the  size  or  shape  of  the  system,  and  by 

the  total  number  of  particles,   N   (unless  otherwise  specified, 

we  deal  with  a  one-component  fluid).   In  actual  systems,  there 

are  always  inhomopeneities  due  to  boundaries  and  to  gravitational 

body  forces.   For  sufficiently  large  systems  however  the  boundary 

inhomogeneities,  being  confined  to  a  region  very  close  to  the 

surface,  i.e.  involving  distances  comparable  to  the  range  of 

intermolecular  forces,  may  be  neglected  when  considering  bulk 

properties  of  the  system  (or  may  be  investigated  explicitly 

by  changing  the  boundaries).   The  inhomogeneities  due  to 

gravity  (or  other  body  forces)  are  usually  very  small  over 

molecular  distances  and  are  treated  by  considering  the  fluid 

as  made  up  of  homogeneous  parts  with  differing  intensive 

parameters.   One  of  the  purposes  of  the  present  paper  is  to 

investigate  and  find  from  the  point  of  view  of  statistical 

^Supported  by  the  AEC  Computing  and  Applied  Mathematics  Center, 
Institute  of  Mathematical  Sciences,  New  York  University,  under 
Contract  AT( 30-1 ) -llj.80  with  the  U.S.  Atomic  Energy  Commission, 
and  by  the  United  States  Air  Force  under  Contract  No. 
AF-I4.9 ( 638 ) -753  monitored  by  the  Air  Force  Office  of  Scientific 
Research  of  the  Air  Research  and  Development  Command. 
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mechanics  a  justification  of  this  procedure.   We  consider 
both  the  case  of  density  variations  alone,  where  we  may 
have  a  true  equilibrium  system,  and  the  case  of  temperature 
and  velocity  variations  as  well,  where  the  system  is  to 
zero  order  in  a  state  of  local  equilibrium. 

In  order  to  carry  out  this  investigation,  we  first 
establish  some  results  concerning  the  form  of  the  joint 
distribution.   n   .(r1,...,r  +.)      of  q+£  particles  in 
equilibrium,  when  the  set  of   q  particles  is  "very  far" 
from  the  set  of  Ji       particles.   The  distribution  nQ+o 
approaches  the  product  of  the  distributions   n   and  n^ 
plus  a  correction  term.   It  is  the  form,  Eq.  (2.22),  of 
this  correction  term  which  we  find  here  under  certain 
conditions.   It  turns  out  that  this  correction  is  related 
to  the  integral  over  all  space  for  the  Ursell  functions  in 
an  infinite  system.   In  the  course  of  the  derivation,  we 
also  find  the  general   l/N  dependence  of  the  low  order 
distributions  for  a  completely  uniform  system  i.e.  one 
confined  to  a  box  with  periodic  boundary  conditions. 
This  case  was  also  studied  by  Oppenheim  and  Mazur  [1]  who, 
"using  a  virial  expansion  and  making  essentially  the  same 
physical  assumptions  as  we  do,  find  a  density  expansion  for 
the  coefficients  of  the  powers  of   (l/N)   in  the  distribu- 
tions  of  a  uniform  system.   It  is  possible  to  show  that 
their  density  expansion  may  be  summed  to  yield  our  result. 
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The  general  correction  terras  to  the  pair  distribution 
function  n2(r,,r2)   when   |r,-r2l  -*-  oo   have  previously 
been  investigated.   When  both  r,   and  r-  are  in  the 
interior  of  a  fluid,  the  asymptotic  form  of  n~   is 
related  indirectly  (through  normalization,  see  e.g. 
Appendix  B)  to  the  scattering  cross-section  of  visible 
light  (of  wavelength  long  compared  to  the  length  of 
molecular  correlations)  by  fluids  found  by  Ornstein  and 
Zernike  [2].   When  the  dependence  of  the  one-particle 
density  n(r)   on  r   is  neglected,  and  some  further 
assumptions  are  made  (implicitly)  which  will  later  become 
clear,  then  [2]  through  terms  of  relative  order   l/N, 

n2(r12)  3  n2g(r12) >  n2[l-nkT#N]  .     (1.1) 

r12-*oo 

Here   n  is  the  average  density  N/V,  X      is  the 
isothermal  compressibility 


y  _   1  Jn 
A  ~  n  Jp 


(1.2) 
T 


and  p   is  the  pressure.   Recently,  one  of  us  [3]  also 
investigated  by  a  very  different  method  (based  on  a  new 
virial  theorem  for  total  momentum  fluctuations)  the 
rather  special  case  of  the  asymptotic  form  of  n0(r, »r0) 
when  r,   is  adjacent  to  a  rigid  wall  which  forms  the 
boundary  of  the  fluid  and  r2   is  far  in"  the  interior. 
It  was  found  there  that,  again  through  order   l/N, 
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n2(^1^2)  r12-»co»   nw  n(1"  I  ^    >  (1*3) 

where 

n   =  p/kT  =  n(r,  )  (l.k-) 

was  shown  to  be  the  value  of  the  density  at  a  rigid  wall. 
Eqs.  (1.1)  and  (1.3)  turn  out  to  be  special  cases  of 
Eq.  (2.22)  derived  in  this  paper.   Another  consequence  of 
Eq.  (2.22)  is  the  vanishing  of  l/N   correction  terms  to 
all  distributions  at   T  =  0.   This  was  conjectured  by 
Peynman  and  Cohen  [ljj.   It  is  here  seen  to  rest  upon  our 
fundamental  hypothesis  of  finite  correlation  length. 

The  basic  assumption  which  is  made  implicitly  in  the 
derivation  of  (l.l)  and  in  Ref.  [1],  and  was  explicitly 
made  in  Ref.  t3J,  may  be  stated  loosely  as  the  absence 
of  long  range  correlations  in  a  fluid.   A  rather 
stringent  formulation  of  this  assumption,  which  in  one 
form  or  another  is  universally  accepted  in  the  theory 
of  fluids,  states  [5,6]  that  if  the  fluid  is  disturbed 
somehow  in  a  limited  region  of  space  of  volume  oj   , 
then  properties  of  parts  of  the  fluid  sufficiently  far 
away  from  u)      (the  effective  unit  of  length  being  the 
range  of  molecular  forces)  will  be  affected  only  to  the 
extent  of   0(u>/V).   Further,  the  effect  of  the  disturbance 
is  supposed  to  approach  its  asymptotic  value  exponentially 
fast.   (The  additional  assumption  of  exponential  approach 
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is  not  necessary  for  many  applications,.   It  is  usually 
sufficient  that  the  approach  go  as  some  power  of  the  inverse 
distance  from  the  location  of  the  disturbance .   What  does 
appear  necessary  to  assume  is  the  absence  of  fluctuating 
behavior  extending  over  the  whole  container. )   It  is  clear 
that  the  assumption  of  no  long  range  correlation  can  be  made 
precise  only  in  the  limit  of  the  volume   V  becoming  infinite, 
the  density  remaining  constant.   It  is  also  clear  unfortunately 
that  a  proof  of  this  central  hypothesis  would  be  extremely 
difficult.   Since  this  hypothesis  is  believed  violated  in 
the  solid  state,  any  proof  of  it  would  give  a  criterion 
for  phase  transitions.   There  is  however  an  ample  experi- 
mental and  intuitive  basis  for  this  characterization  of  the 
fluid  state^  we  shall  accept  it  here  and  use  it  in  our  proofs 
in  this  paper,  leaving  a  more  formal  discussion  to  Appendix  A. 
We  give  two  derivations  of  our  general  result,  a  thermo- 
dynamic one  is  in  Section  2,  and  a  statistical  one  based  in 
part  on  the  virial  expansion,  and  thus  presumably  valid  only 
for  gases,  in  Section  3»   The  application  of  our  result  to 
non-uniform  systems  is  presented  in  Section  L|_„   In  Appendix  B, 
we  extend  the  result  on  asymptotic  expressions  to  mixtures, 
and  also  consider  light  scattering  from  mixtures.   Appendix  C 
contains  a  formulation  of  our  main  result  in  terms  of  Ursell 
functions,  and  Appendix  D  contains  a  further  self -consistent 
confirmation  of  this  result. 
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2.   Thermodynamic  Derivation. 

In  order  to  make  our  arguments  concrete,  we  shall  consider 
explicitly  the  asymptotic  value  of  the  pair  density  n^r-^r,,). 
Our  argument  however  will  be  of  such  a  form  that  extension  of 
the  results  to  the  other  distributions  will  be  immediate. 
From  its  definition,   n2   is  symmetric  in  r,   and  j?2.   When 
np(r-i,r.?)   is  integrated  over  the  whole  pair  space,  each  pair 
is  counted  twice,  thus  yielding  the  normalization 

J  [  n^r^r^)  dr^d^  =  N(N-l)  (2.1) 

where      N     is   the   number   of   particles   and  the    integration  is 

over   the  system   volume      V.       (Since      n2(r,,£2)      is   zero  when- 
ever    r,      or     Y~      is   outside      V,      the    integration  can  be 

extended  over   all   space,   which   is   indeed  necessary  when  there 
is   no  well-defined   volume      V. ) 

The  assumption  of  a  finite   correlation  length   implies, 
for     r,~     very   large, 

n2(rltv2)      >     n(r1)    n(r£)    [1   -|(l+o(l))]    .    (2.2) 

To  show  this,  we  first  write  Eq.  (2.2)  in  a  form  which 
exhibits  more  transparently  its  relation  to  the  correlation 
length.   Consider  the  conditional  number  density 

w(r2|r1)   =  n2(r2,r1)/n(r1)         .(2.3) 

for  particles  at  r~  when  it  is  known  that  there  is  a 
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particle  at   r, «   Eq„  (2„2)  may  then  be  written  as 

w(p0|rn)   — ►  n(p«)  [1  -  b/N  +  .00]  .   i  .  '(2.1|.) 

Now  in  a  classical  system  in  equilibrium,  the  presence 
of  a  particle  at  r,   does  two  things,  (l)  it  leaves  only 
N-l  particles-  with  unspecified  positions  and  (2)  it  intro- 
duces an  effective  one-partiole  potential  fi^M^     for  ^e 
j**1  particle  (here  0(r)      is  the  assumed  potential  between 
any  two  moleculesh   The  conditional  density  w^glr^)   is 
therefore  identical  with  the  ordinary  density  "in  a  system 
of  N-l  particles  with  this  extra  perturbing  potential.. 
Hence  for  distances  large  compared  to  X,   the  range  of 
0s     w^2^1^   should  be  equal  to  the  density  at  r2   in  a 
fluid  consisting  of  N-l   particles,  plus  a  term  O(A-yv). 
Since  the  difference  in  the  density  at  r^     between  a  fluid 
consisting  of  N  particles  and  one  of  N-l   particles  is 
again  0(1/N),   we  are  led  to  (2„ljJ.   For  the  special  case 
of  an  ideal  gas,  $   =  0,   and  we  find  immediately  that 

w^2'Sl)   =   (N-VN)  n(r2)   =  n(r2)[l-  l/Nj.    (2.5) 
Prom  its  definition,  we  of  course  always  have 

fwd^lr-j.)  ^2   =  N"10  (?"°6) 

For  a  real  fluid,  the  added  "external"  potential 
0(r„,  )   produces  a  change  in  density  in  the  vicinity  of  ^r, 
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with  the  result  that  the  effective  number  of  particles 
whose  position  is  unspecified  (which  determines  w2  far 
away  from  r,  )   differs  from  N-l,   and  (2.5)  is  replaced 
by  ( 2 . It.) .   We  can  express  this  idea  more  clearly  in  the 
following  way.   Let  us  divide  the  volume   V  into  two 
macroscopic  parts   VA  and   Vg,   VA  +  VB  =  v*   The  1Part 
VA   contains  the  point  £,  ,   while  all  points  r^     in  VB 
are  sufficiently  far  from  r,   that   n2  may  be  given  its 
asymptotic  value.   The  average  number  of  particles  normally 
in  V.   and  Vg  will  be  designated  by  N"A  and  Ng,  where 
NA   and  N_   are  regarded  here  as  both  of  order  N.   Now 
for  any  particular  partition  of  the   N  particles  into  NA 
in  VA  and  N„  in  Vg,   where 

Ng  =  N  -  NA  ,  . (2.7) 

we  may  write 

n2(£l,r2|NA,NB)   =  n(r1 |r2,NA,Ng)  n(r2|NA,Ng).     (2.8) 

Here   n(r,  l£p>NA>NR^   *s  ^ne  densit7  at  .J-i   when  there  are 
NA  particles  in  V.   and  NR  particles  in   Vg,   one  of  the 
latter  being  specified  to  be  at  ^2;   similarly, 
n(r2|NA,Ng)   is  the  density  at  r?  when  there  are  N. 
particles  in  VA  and  NR  in  Vg.   This  situation  corres- 
ponds to  a  model  system  in  which  V.   is  separated  from  Vg 
by  a  surface  S.     which  is  impenetrable  for  particles  cross- 
ing between  VA  and  Vg,   but  which  does  not  otherwise 
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modify  the  Hamiltonian  of  the  system . 
We  now  have  generally  that 

n2(r19r2)     =  HN  P(HA)  n^T^Tl-S^,  (2.9) 

where     P(N.)      is    the   probability  of  having     NA     particles 
in     V.,      and  correspondingly     N-NA      in     Vg.      Thus  performing 
a  Taylor  expansion  of     n2     about      NA  =  N^, 

+  |<(NA=NA)2^(a^Nf)[n(r1lr2,IA,N-NA)n(r2|NAPN-NA)]      (2.10] 


"•       o  o  o      5 


where      <(NA-ITA)2>    =     HNfl   P(NA)    (VHA)2    ' 

We   assert   that    (l)    the  Taylor  expansion  may  be    stopped  at 
the   term     <(NA~NA)2>     and      (2)      nCrJjr^S^N-l^)      in   (2.10) 
may  be   replaced  by     n^lN^N-N^).      The    first  point   follows 
from  the   fact  that,    judging  from  the   prototype    of   an  ideal 

gass     <INA~^aO    =     0(^A    ^      Slnce   the   volume      VA      is 
macroscopic,    the  remaining  terms   in  the    series   are   of  order 

|pV2   q(i/na)      and  may  be    dropped.      Secondly,    the   fixing  of 

A  A 

a  particle  at  r-   in  Vg  for  given  values  of  NA  and  Ng 
can  affect  the  density  at  r,   in  V"A  only  by  modifying  the 
distribution  of  particles  in  Vg  close  to  the  surface  2>A'9 
this  in  turn  modifies  the  neighboring  distribution  in  V"A 
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close   to     S.,      which  can  then  affect   the  density  at    £-1  • 
Since   the  potential   change   due   to    the   positional  perturba- 
tions  of  the   particles   near  the    surface   in     VR     extends 
only  to  the   order   of     X      inside  region     V.,      the    assumption 
of   short-range   correlation  tells   us   that  when    j?,      is  far 
from  the    surface      SA, 

nfr-Jr^N^N-iy      =     ntrJI^N-l^)    + 

(7lSa/Va)  [order   of  deviation  of  n(rg  |r2,N"A,N-NA)   (2.11) 

from  n(rs|NA,N-NA)J  , 

where     jr.,      denotes    a   surface      particle    in     VR.      But    again  by 

the   assumption  of    short-range    correlation,    the    deviation  in 

(2.11),    elicited  by   a  perturbation  over  a  volume     ~-»  X   ,   will 

be   of   order     0(X3/Vg)    ^  0(l/N)   =    (NA/N)0(l/NA) .      The 

correction  term  in  (2.11)  for  r,   and  r0  far  from  S. 

—  1  ~2  A 

will  then  be  /^-»n(XSA/NA)0(l/N)      and  will   vanish  compared  to 
0(l/N)      as      V.      increases. 

The  foregoing   argument    can  be    improved  by  further 

t  " 

dividing   V„  into   V"B,   surrounding  r-,   and  V*B   separat- 

t 

ing     VR     from     V"A.      We   shall  however  not   attempt   this   here. 

Accepting  then  the    above   argument,   we   now  rewrite   Eq.    (2.10) 
asymptotically  as 

n2{^1^2)      -*      n(£llfA'N-V    n(i:2'IA'N-\) 

(2.10') 

+  |(NA-NA)2    {d2m2k)    [n^lN^N-^)   n(r2|NA,N-NA)]. 
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The  steps  which  led  to  (2.10)  may  also  be  utilized  to  write 
n^)   =  EPtN^^ly-IJ  (2.12) 

=  ntrjN^N-^)  +  f^V^)2^2/^?2)^  |NA,N-NA) 
n(£2)   =  n(r21!A*N"WA)  +  |((NA-fA)2/>(^2/^f)n(r2|NA,N-NA). 

In   (2.12),   we   have   neglected  terms   of   smaller   order  than 
1/N.,      but    this    is    consistent  with      V.      being   of  macro- 
scopic   size.      Combining  Eqs.    (2.11)    and    (2.12)    then  yields 
in  the   asymptotic  region 

n2(£l>£2}      ~ *     n(£l}    n(?2)  (2.13) 

+  <(6NA)2^/dNA)n(i;1|NA,N-NA)(^/^NA)n(r2|NA,N-NA) 

where      6N„    =  NA-Nft . 
A  A     A 

Since  the  second  term  in  (2.13)  is  already  of  order 
l/N,   any  parts  of  it  which  are  of  further  order   o(l) 
can  be  omitted.   This  permits  (since   p-   and  v~     are  far 
from  the  surface)  the  replacement  to  order   S/V 

ntrJI^N-^)  =  nfe-JI^o)  ^^ 

n(£2|NA,N-NA)   =  n(r2|o,N-IA)  . 

Both  N.   and  NR  are  functions  of   N.   We  may  therefore 
express  (2.13)  in  terms  of  derivatives  with  respect  to   N. 
Thus 
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^n(r2|o,N-NA)  dn(r2|o,NB) 


C^N,  dN 


B 


dN  x    dN' 

^n(^2}    dN 
"'T*N-B 

and   similarly 

^(rjN^o)  ^n(gl}    dN 


dn(rJo,Nn(N))      dIL     , 

-,   (-J)"1  (2.15a) 


(2.15b) 


again  neglecting  terms    of  order      S/V. 

Setting     n  =  N/V     and   leaving  the   n-dependence    implicit, 

Eq.  (2.13)  thus  reduces  to 

,        ^n(r,  )       e)n(rj  ? . 

.ifdlL+M.]  (2-16) 

dNA        dNB 

where  we  have  used  the  identity   ( dN/dN . ) ( dN/dNg )  = 
dN/dfA  +  dN/dNB. 

We  now  show  that  the  expression  (2.16)  is  independent 
of  the  division  into  volumes  V.   and   V"B,   and  find  its 
explicit  form.   It  follows  from  the  general  principles  of 
statistical  mechanics  that 
P(NA)  =  exp-[FA(T,VA,NA)+FB(T,VB,NB)-F(T,V,N)J/kT 

C     !  (6NA)2   ^FA(T,VA,ffA)   a2FB(T,VB,NB)^>^) 
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F  being  the  Helmholtz  free  energy,  and  we  have  again 
consistently  neglected  terms  of  order   o(l)„   Here   N.   is 
determined  from  the  condition 

<)F A(T,VA,flA)  ^PB(T,V  N-N  ) 

n  -  — ^—  =  ^b  =  9« =  > . 

and  |X  =  nA(NA(N))  =  (i(N)   is  the  chemical  potential  of  the 
full  system.   This  readily  yields  the  expression 


<(6NA)2>  o  kT^  (<^_  +  diL)]"l  (2.18) 

or,    combining  with  Eq»    (2.16), 

kT         Snir-,  ) 

n2(ii^2)    -*   n(^i)n(^2)  "  T  [n  —&r]  * 

,m^£i  r1^,  (2'19) 

For  a  uniform  system  with  no  body  forces,  the  quantity 
(l/N)^N/9p,  is  equal  to  n  times  the  isothermal  compressi- 
bility /-<,   For  a  general  non-uniform  system,  we  may  still 
write 


where  ^  is  now  some  average  compressibility <»   This  leads 
us  finally  to  the  asymptotic  expression 

n_(r.,r5)    ■»  nfr.WrJ-^ninfrJJlninfrJ].    (2.21) 
2  ~\P'-<2  ~1        *-2  N  ^n        "-1  an       ^2 

When     r,      and     r0      are   both   in  the   interior   of   a 
-^1  ~2 

uniform  fluid,    (2.21)   becomes   the    Ornstein-Zernike  relation 
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g(r)      — *     1  -   nkT^/N   ,  (1.1«  ) 

while    if     r,       is   at   a  rigid  wall  then  as    shown  in  Ref.    [3 J, 

n       =     n(r.    at  wall)      =     p/kT, 

cMr-^/an     =     V(nkT/)    , 
and 

n2(£l^2}      -*      nwn[1   "   n^lJ    »  (l-3') 

a  result   derived   in  Ref.    [3J. 

The   extension  of  the   above   analysis    to   the   distribution 

of     q+j£    particles   is    straightforward  and  yields 

<)nn  dnt 

n 


W 


»,r¥^'^i. 
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3«   Statistical  Mechanics  Derivation. 

We  shall  now  investigate  the  asymptotic  behavior  of 
n  .„  from  a  statistical  mechanical  viewpoint.   To  do  this 
involves  two  steps;   (l)  a  decomposition  of  quite  general 
validity  of  thermodynamic  quantities  and  distribution 
functions  into  a  dominant  part  and  an  0(l/N)   correction 
term  expressed  in  terms  of  the  dominant  part.   (2)  a 
computation  of  the  spatial  asymptotic  form  of  the  dominant 
part  of  a  distribution  and  consequently  of  the   0(l/N)   term, 
in  which  the  explicit  form  of  the  canonical  partition 
function  plays  a  crucial  role.   Since  the  dominant  part  in 
fact  represents  a  grand  canonical  average  from  which  the 
fixed-N  components  are  then  extracted,  there  is  a  family 
relation  with  the  "box  within  a  box"  thermodynamic  argument 
of  the  previous  section. 

The  equilibrium  statistical  mechanics  of  an  interacting 
N-particle  system  in  a  volume   V  may  be  obtained  from  the 
partition  function  QN,   or  from  its  physical  counterpart, 
the  free  energy 

PN  =   -CL/p)  In  QN  ,  (3.1) 

where  kp   is  the  uniform  reciprocal  temperature.   Q„   is 
given  explicitly  by 

^(1/N1)  J"  exp  (-pHjjHdr  dp/h3)N  ,   classical 

Tr  (exp  -0%)  ,  quantum  mechanical, 
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Hw     being   the    system  Hamiltonian.      However,    the   precise   form 
of     QN     is   unimportant      at   this    stage.      What   is   vital   is  the 
smallness   of   density  fluctuations,    to    the   extent    that   the 
grand  partition  function 

Q(z)     =     2Zo°°   QM  zM  (3.3) 

at  fixed   V  and   z   has  contributions  from  only  a  very 
•small  range  of  particle  number  M.   If  this  is  so,  we  may 
contemplate  a  Darwin-Fowler  or  steepest  descent  extraction 
of  Qjj  from  Q(z). 
Let  us  define 

L(z)   =  In  Q(z)  .  (3.i4.) 

For  an  open  system  or  grand  canonical  ensemble,   L   is 
related  to  the  pressure  by  L  =  pV/kT  and  z   to  the 
chemical  potential  by   P(x  =  lnz.   Here,  the  quantity  L 
is  but  a  means  to  an  end,  for  we  may  express  QN  in  terms 
of  it  through  the  contour  integral 

QN  =   (l/27ti)  ^exp  [LU)J  £"N  &Z/Z, 

(3.5) 

=      (l/27t)   f  exp  [L(re19)-Ni9-N  In  rj  d9  . 
J-x 

The  large   N   asymptotic  form  of  (3.5)  appears  from  a 

steepest  descent  evaluation  as-^ 

1 
%     2?  z"N  eL(z)  [2a  zjL  z  jLL(z)]~2     (3#6) 
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where 

z  ^L(z)/^z     =     N  (3.7) 

determines      z.      If  the  method   is  valid,    corrections   to    (3-6) 
for   derived  quantities    such  as     F„     may  be    shown  to   be   of 
order     0(l/N    )      compared  to   the  principal   term,    and   this   is 
strictly  negligible    for  most   applications,    including  ours. 
Such  a   statement   is   as   always   fully  meaningful   only  for   a 
uniform  fluid,   where   one   can  take   the    limit     N   —*■   oo       at 
constant   density,    but   the    same   numerical   order   of  magnitude 
should  hold  for   a  non-uniform  system.      Prom    (3.1),    we   have 
as  well 

-  pFN     =     L(z)    -  N   In  z   -  -|  In  2* 

(3.8) 

On  the  basis  of  (3*8),  we  may  introduce  the  "dominant 
part"  of  FN, 

-  pF^0)   =  L(z)  -  N  In  z  -  |  In  2%   .  (3.9) 

f!  s    which  coincides  with  the  result  of  the  usual 
Mayer-Ursell  virial  expansion  [8],  is  of  course  just  the 
grand  canonical  average,  and  would  yield  the  correct  free 
energy  per  particle  in  an  infinitely  large  uniform  system. 
In  a  general  system,  the  principal  property  of  FN 
which  will  be  utilized  is  the  following.   Consider  any 
variation  of  the  system  which  does  not  affect  total  particle 
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number.   Then  according  to  (3»9), 


N 
or  inserting  (3»7), 


-  6(p  P^0))   =   6L(z)  -  (N/z)  6z  , 


-  6(p  F^0))   =  5L(z)  -  [^L(z)/^z]  6z  . 
Hence  we  may  write 

-5(PP<0))   =   6L(Z)lz=const>         (3.10) 

Eq.  (3.10)  permits  us  to  compute  a  variation  of  F„ 
explicitly  in  terms  of  L  without  being  required  to  solve 
for   z.   Finally,   F„   is  to  be  expressed  in  terms  of  F 
Since 

-  p  <}F^0)/9N  =   (^z/^N)(^L/(3z-N/z)  -  In  z  =   -In  z, 
i.e.  the  thermodynamic  relation  |x  =  <3F/^N,   then 
P  ^2F^0)/c^N2  =  (l/z^z/^N  =  (z  ^N/^z)"1  =  (z  <^z  z^z  L)"1. 

Inserting  into  Eqs.  (3.8)  and  (3«9),  we  have 

d2F{0) 
PFN  =  Pp*0)  -|m  tp-^5-).        (3.11) 

Eq.  (3.11)  may  be  rewritten  in  terms  of  the  compressi- 
bility.  For  a  very  large  uniform  system,   F^/N  is  a  function 
of  n  =  N/V  alone.   Hence  the  isothermal  compressibility  may 
be  written  as 
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X  =    (n  Jv/dn)"1  =    (-N  d/dK  ^F/av)"1  =    (-N  ^N  N  ^/^V  P/N)"1 

=    (n  d/d®  N2  ^/^N   P/N)"1   =    (n  N  ^F/^N2)'1    . 

For   an  arbitrary  system,   we  may   correspondingly  define   an 
average    compressibility 

X  -    <|»?-4fJ0,>-1  <3.12) 

dn 

where  differentiation  is  at  constant  volume  (or  constant 
external  potential) .   Eq.  (3«12)  is  seen  to  be  equivalent 
to  the  expression  (2.20).   Further,  applying  a  variation 
6,   commuting  with  N,   to  (3. 11)  yields  by  virtue  of  (3.12) 

__      _  (0)      1  nX  -2    J*  B  P(0)  M  n) 

our   fundamental  relation.      The   second  term  on  the   right 
hand   side   of    (3.13)    is   a   correction,    containing  the    factor 
l/N.      For   a  large   uniform   system,    the   relative    correction 
(FN      itself  goes    as      N)      is   literally  a     l/N      term,   while 
for   a   non-uniform   system,    this   is    still   a  reasonable    assess- 
ment  of   its  numerical   value. 

Relation   (3.13)    is   readily  applied  to   determining  the 
N-dependence   and   asymptotic  form  of  the    configuration  space 
distribution  functions.      For   the   N-dependence,    it   suffices 
to  express   the   distributions   in  terms   of  the   free   energy. 
Suppose   that   the    system  Hamiltonian  contains   an  s-body 
potential, 

H  =   ...   +   (1/si)  2Zj.  j£ooe/±     &  (r-i»00»£i    )+"°»        (3.U4.) 
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which  in  the   final  evaluation  is   to  be   given  its   actual 
value,    perhaps   zero,    but  not  until   then.      Then   it   is   easily 
shown  that   the    s-body  distribution  function     ng,      normalized 
so   that 

fn   (r,,...,rj    dr,,...,dr        =     Nl/(N-s)l    ,  (3.15) 

v        S  /»\J.  VS  ~X  ~vS 


may  be  written  as 


V^i"-"^ 


si  JL 

P  ft, 


iQN 


N  60 


(77 


(r,, . . .,r  ) 

~>-J.      "•  s 


or 


r(s) 


ns  =   si  6FN/6^vs;  .  (3.16) 

By  making  0  other  than  a  pure  configuration  space 

potential,  this  result  extends  at  once  to  phase  space  distri- 
butions and  density  matrices.   Now  applying  (3.13)  to  (3.16) 
results   in  the  relation 

,2 


n 


(0)    1  nX  n2  Jt      (0) 


(3.17) 


where   n   ,   arising  via  (3.16)  from  FN   ,   actually 
represents  a  grand  canonical  average  distribution.   When 


n 


(0) 


depends   only  on     n  =  N/V,      then  it   coincides  with 


n       in  an  infinite    system,    and    (3*17)    expresses   the      l/N 
corrections.      For   example,   with      s  =  2      and  large    separation 
r,p,      then  in  a  uniform   system,    (3.17)   reduces   to  the 
previously  discussed     n~  =  n   (l-n%/pN).      For   a  non-uniform 
system,    the   utility  of    (3.17)    derives   from  the   special 
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properties  which  n     possesses,  and  in  particular  its 

5 

simple  asymptotic  form. 

Our  major  objective  is  to  find  the  asymptotic  form 
of  n  +.  when  the  coordinates  separate  into  two  distinct 

bunches 

q  =  (^i'---^q)        and       -^  =  (jrq+i"-"iW  • 

For  this  purpose,  the  corresponding  asymptotic  form  of 
n  +L      is  required,  and  it  may  be  obtained  by  assuming 
sufficiently  rapid  convergence  of  the  standard  cluster 
expansion  for   L(z),   a  stronger  condition  than  the 
previously  employed  finite  correlation  length.   The  point 
is  this.   First   consider  classical  equilibrium  statistical 
mechanics.   If  the  factor 

f(s)(r,,...,r  )   =  exp  (-p0( s } (r, , . . . ,r  ) ) -1      (3.18) 

is  regarded   as   a   "star"    connecting  vertices     r,,...,r    ,      and 
if  we    further    introduce    the   Ursell   factor 

\(rlt...tvt)     =     SlDt    irls...tvt)  (3.19) 

where  the   D.   run  over  all  (non-repeated)  products  of  f 
represented  by  diagrams  in  which  the  r~  ,  .  .  .  ,r .   are 
connected,  then  it  is  known  (see  Hill  [2],  p.  136)  that 
Q(z)   of  (3.2)  and  (3-3)  achieves  the  form 

Q(z)   =  exp  ^2   t  f Ut(r1,...,rt)  d*r  zt   Gt/tl]    .     (3.20) 
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Here   G   is  a  factor  resulting  from  momentum  integrations. 
Prom  (3.16)  and  (3.10),  then 

n(0)   =  £-  zt  ct  (.sl/p)  (5/6)ZJ(s))J)Ut  d^/ti  ,   (3.21) 

whose  simplicity  is  due  to  the  fact  that  explicit 
z-dependence  is  untouched  by  the  variation. 

Suppose  now  that  (3.21)  converges  rapidly  enough 
that  only  terms  up  to  some  finite   t   need  be  considered. 
For  finite  range  forces,  this  then  requires  only  clusters 
up  to  some  finite  spatial  diameter,  related  of  course  to 
the  correlation  length  but  conceivably  depending  on  s. 
Sets   q  and  ^/   in   s  =  q+J     are  then  to  be  regarded  as 
asymptotically  separated  if   q  and  JC,    are  separated, 
member  from  member,  by  more  than  this  finite  diameter. 
If  under  these  conditions   6jDt  d*r/C0*q   '   did  not 
disconnect  set   q  from  set  Jl       in  D.  ,   the  cluster 
would  have  an  impermissibly  large  diameter.   Thus  a 
separation  must  occur.   We  note  that  U./tl   lists  all 
distinct  diagrams  (with  undesignated  vertices,  but 
weighted  by  the  number  of  repetitions  under  vertex  permuta- 
tion), while   (-sl/p)6/6jZJls^   for  separated  q  and  JL 
removes  a  "star"   f^3'   in  any  order  in  which  it  may  be 


presented.   It  follows  that,  including  both  possibilities 
f  q   and  1  for  the   q  conned 
the  JL ,   the  separation  is  into 


f  q   and  1  for  the   q  connection,  and  similarly  for 
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properties  which  n     possesses,  and  in  particular  its 
simple  asymptotic  form. 

Our  major  objective  is  to  find  the  asymptotic  form 
of  n  ..   when  the  coordinates  separate  into  two  distinct 
bunches 

9  =  {r_l>'->Zq]  and      ^  =   (JTq+l"--^qV)  * 

For  this  purpose,  the  corresponding  asymptotic  form  of 
n  +i      is  required,  and  it  may  be  obtained  by  assuming 
sufficiently  rapid  convergence  of  the  standard  cluster 
expansion  for   L(z),   a  stronger  condition  than  the 
previously  employed  finite  correlation  length.   The  point 
is  this.   First   consider  classical  equilibrium  statistical 
mechanics.   If  the  factor 

f(s)(r,,...,r  )   =  exp  (-p0( s ] (r, , . . . ,r  ) ) -1      (3.18) 

is  regarded   as   a   "star"    connecting  vertices     r -.,..., r    ,      and 
if  we   further    introduce   the  Ursell   factor 

Vsi,.v.fpt)   =  21  Dt  ^i'--'^  (3-19) 

where  the   D.   run  over  all  (non-repeated)  products  of  f 
represented  by  diagrams  in  which  the  r~,...,r.       are 
connected,  then  it  is  known  (see  Hill  [2],  p.  136)  that 
Q(z)   of  (3.2)  and  (3>3)   achieves  the  form 

Q(z)   =   exp  21  C  fUt(£l'"''£t^  ^  Z%   ct/tJJ  •  (3.20) 
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Here      C      is   a  factor  resulting  from  momentum  integrations. 
Prom    (3.16)    and    (3.10),    then 


JO)      _     Y~     t 
ns  -     2_  z 


C*    (-sl/p)    (6/60(s))Jut   c^r/tl    ,       (3.21) 


whose  simplicity  is  due  to  the  fact  that  explicit 
z-dependence  is  untouched  by  the  variation. 

Suppose  now  that  (3*21)  converges  rapidly  enough 
that  only  terms  up  to  some  finite   t   need  be  considered. 
For  finite  range  forces,  this  then  requires  only  clusters 
up  to  some  finite  spatial  diameter,  related  of  course  to 
the  correlation  length  but  conceivably  depending  on  s. 
Sets   q  and  jl      in   s  =  q+J     are  then  to  be  regarded  as 
asymptotically  separated  if   q  and  &    are  separated, 
member  from  member,  by  more  than  this  finite  diameter. 
If  under  these  conditions   6  J  Dt  6.^/60   q  '      did  not 
disconnect  set   q  from  set  j2   in  D.  ,   the  cluster 
would  have  an  impermissibly  large  diameter.   Thus  a 
separation  must  occur.   We  note  that  U./tJ   lists  all 
distinct  diagrams  (with  undesignated  vertices,  but 
weighted  by  the  number  of  repetitions  under  vertex  permuta- 
tion), while   (-s!/p)6/60(s^   for  separated  q  and  JL 
removes  a  "star"   f     in  any  order  in  which  it  may  be 


presented.   It  follows  that,  including  both  possibilities 
f  q   and  1  for  the   q  conned 
the  JL ,   the  separation  is  into 


f  q   and  1  for  the   q  connection,  and  similarly  for 
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(-qj/p)(6/60(q))f  Dm  drnr/mj(-/£i/p)(6/6^U))JDndnr/nI  . 

Here  ra+n  =  t,   and  all  distinct  combinations  D  ,  D   are 
'  m*   n 

obtained  in  this  fashion.   Inserting  into  (3.21),  we  at 
once  conclude  that 

n(0)   _^    (0)   (0)  ,-  2-j 

for    asymptotic   separation  of   sets      q      and     j2  .      This  may 
be  regarded  as   the    characteristic  property    of     n        : 
the   grand   canonical    distributions,    for    asymptotic 
decomposition  into   particle    subsets,    separate   through 
order     l/N     into   their    component   distributions.       [It  might 
appear   that   a  corresponding   statement   could  be  made   for   the 
diagrammatic   density   expansion  of   a   canonical   ensemble  with 
periodic  boundary    conditions   and  translation-invariant 
potential,    but    since   the   generating  function  for   the 
irreducible    clusters    is   no   longer  the    free   energy,    and 
since   more   diagrams  may  be    decomposed  by  the    variational 
differentiation,    the    form    (3*2.2)    is   not   obtained.      Instead, 
one    finds    (3«23)    directly,   but  more  laboriously.] 

Applying  the    correction  formula    (3.17)    to    the    combina- 
tion    n    ..  -  n     ru    ,      we  have 

-   [n(0)-(n^2pN)d2n'0)/an2J[n<0)-(n3f/2PN)d2n(0,/dn2J    . 
q  q  a> 
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p 

Employing    (3.22),    and  dropping  the   final      0(l/N    )      term   as 
well   as   the   deviation   of     n        from     n  within  a   correction 

term,    then 


n 


q+z  —  «\,  it  -  *  n£  t»  &'  (n&-        (3-23) 


the  desired  asymptotic  relation.   The  foregoing  is  classical. 
For  quantum  mechanical  distributions,  the  explicit  form  of 
the  Ursell  factor  is  altered  [9],  but  similar  comments  are 
appropriate.   The  non-commutativity  of  coordinates  and 
momenta  however  introduces  additional  effective  coupling 

in  the  form  of  propagation  factors  of  range  of  the  order 

P    1  /p 
of  the  thermal  de  Broglie  wave  length,   ((3h  /2m)  '  , 

thereby  contributing  to  the  maximum  effective  cluster 

diameter.   At  very  low  temperature,   p  — >   oo ,   greater 

care  is  required,  as  it  is  when  infinite  range,  e.g. 

Coulomb  forces,  are  present. 

We  remark  here  further  on  the  relation  of  Eq.  (3.17) 

to  the  general  N-dependence  of  the  low  order  distribution 

functions.   This  may  be  of  some  direct  relevance  in  dealing 

with  systems  in  which  N  is  actually  a  small  number,  such 

as  those  used  in  machine  calculations  by  Alder  [10]  and  others. 

When  our  assumption  concerning  the  convergence  of  the  virial 

expansion  is  valid,  then  for  a  uniform  system,   n     depends 

only  upon  the  density,  in  which  case  as  indicated  following 

(3.17),   n     coincides  with  the  leading  term  in  a  l/N 
s 
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development  of  n  .   Eq.  (3»17)  then  states  a  verifiable 

s 

relation  between  the  N- independent  term  and  the  first 
correction  to  it,  which  is  proportional  to   l/N» 

The  relation  (3°23)  may  also  be  extended  to  the 
intrinsic  correlations  or  Ursell  distribution  functions. 
This  is  carried  out  in  Appendix  C„ 
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Ij..   The  Local  Nature  of  Distribution  Functions. 

We  shall  now  apply  the  results  of  the  previous 
sections  to  a  one -component  fluid  whose  intensive 
parameters  such  as  density,  temperature,  and  local  velocity 
vary  slowly  with  position.   The  Gibbs  ensemble  describing 
such  a  system  has  been  considered  by  many  authors  [11]. 
The  classical  ensemble  density  n,(X),   where 
X  =  j . . .r. , . . . ,p. . . .  ?   is  a  point  in  the    -space  of  the 
system,  is  generally  written  as  a  local  equilibrium  part 
\i(X)      and  a  correction  term  [i    (X).   The  part   |x(XJ   is  a 
superposition  of  canonical  (or  grand  canonical)  ensembles 
for  each  small  region  of  the  fluid  at  its  own  temperature 
T(x),   velocity  v(x),   and  density  n(x) : 

[i(X)      =     £(X)  +  H-'U)  (Ij-.l) 

J(X)  =  (l/NlQ)exp-fdx  p(x)[E(X,x)-v(x).p(X,x)  - 

(k-2) 
-   v(x)  n(X,x)J  . 

Here      x      is   a  point    in  physical    space,      P(x)    =    (kT(x))      , 
and     v(x)      is    chosen  to   give   the    correct   density      n(x). 
E(X,x),      P(X,x)      and     n(X,x)      are   the  microscopic   energy, 
momentum,    and  particle   density   at     x     when  the    state   of 
the    system   is   represented  by  the   phase  point     X   [12]: 

E(X,x)   =  Til   SU-r±)[p2±/2m  +  f  Zj   t^ ± ^) +XJ (Sl)  1 

P(X,x)   =  H*   6U-v±)   p±  (k-3) 

n(X,x)   =  T~®  6(x-r.)    . 
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0(r. .)   and  U(r, )   are  the  internal  and  external  potentials 


/\ 


(K'K) 


respectively;   Q  is  a  normalization  constant,  reducing  to 
the  canonical  partition  function  Qw  for  an  equilibrium 
e  ns  emb 1 e .   Thus 

where 
^un.QO  =  exp  -S*{-T^  +  P^i^ 

-[(pi-mv(r1))2/2m+|  JIj  ^(^ij)jj  » 

Q  =   (l/NDJ^^X)  d^  ...  dpN 
and 

r(r±)    =    [-u^J  +  ^mv^  +  v^)]  p^)  .    (k.5) 

The  term  \x  (X)   in  (I4..I)  is  of  first  order  in  the 

gradients  of  the  hydrodynamical  variables  in  the  sense  of 

vanishing  in  the  case  of  uniformity,  and  has  a  vanishing 
integral.   It  is  u.   which  is  supposed  responsible  for  the 

A 

dissipative  behavior,   while      u-      should  yield  the   equilibrium 
form,    at    the    local   value   of  the    intensive   variables,    for   such 
quantities   as   energy   density     E(x),      pressure     p(x),      etc. 
This  will  be   true    if  the   low   order  distributions      fg,      and 
in  particular   the    two-body   distribution   function 
f0(r,  ,r„,p,  ,pJ.      computed  from    £     have   their  equilibrium 
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form  at  the  local  value  of  the  intensive  variables.   We 

shall  investigate  the  validity  of  this  hypothesis. 

Now  the  distribution  f   for  the  position  and 

s 

momentum  of   s   particles  factors  readily  into  a  position 

distribution  n  (r, ,...,r J   and  a  momentum  distribution 
s  ~1*    ~s 

(with  positional  dependence  in  the  large).   We  have 

^s(rp-'is^"'Ps)  =  (Ni/(N-s)1)j£^s+r-%  dSs+r-d?N 


n   (r,,...,r    )   IT? <(2rankT(r . 
s  -1*        *^s  1  /  -^-i 


))3/2<  (I4..6) 


x.  exp  -  [p(r.  )  (p.-mv(r.  )  )  /2m]  (  . 

The  momentum  distribution  is  thus  locally  Maxwellian  and 
we  need  consider  only  n  (r, , . . . ,r  j .   The  position  distri- 
bution  n   is  a  functional  of  the  temperature   T(r)   and 
(through  \)      °£   the  density  n(:r)   in  the  whole  container, 
It  is  our  aim  to  show  that,  for  slow  variation  of   p(r) 
and  n(r),   the   n   depend  only  on  the  values  of  these 
quantities  in  the  region  containing  r,,...,r  .   In  fact, 
it  will  develop  that,  for   s  =  1,2 

A.  

%(£.,_, ...,rs)   =  n3(r1,...,rs;^,n)  (1^.7) 

to  second  order  in  the  gradients.   Here   n   is  the 

s 

equilibrium  distribution  for    a   system   at  uniform 
temperature      (k(3)~        and  density     nj        "J3  =  p(r)      and 
n  =  n(r),      where     r   =    (l/s)  5~f  r,      is   the    s-particle 
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centroid.   For   s  >  2,   gradients  at   r  will  also  be 

required.   However,  there  will  be  no  contributions  from 

outlying  elements  of  the  fluid,  so  that  even  for   s  >  2 

the  distribution  function  is  truly  local,  as  expected. 

We  shall  consider  first  the  change  in  n   due  to  a 

s 

change  in  the  function  y^)-   It  readily  follows  from 
Eqs.  (lj.«I|.)  and  (I4..  6)  that 


(I|..8) 


6n  (rn  ,...r  ) 

S  —1  '    *~  S    _  ^   /  \  TT—  s   c  /         1 

jt- 7 — 1 =  n  (r.,..,r  )   >   6lr,-rJ 

Svlr;        s  ~1*   »— . s  4r-    ~i  ~ 

+  n  .,(r,,..,r  ,r  )-n  (r,  , .  .,r  )n,  (r )  . 
s+1  a.1*  vsv    s  A~3.*      '<~  s     j.  ->^ 

in,  =  n,   is  the  density  which  determines  y)  • 
Hence  under  the  infinitesimal  alteration  6y(r), 

6n  (r, ,..,r  )  =  n  (r,,..,r  )  Yjf  6y(r.  )  +   [n  .,  (rn,..,r .r) 
s  ~^i»   >~.s     s  ~*±'      »^s  * — 1    -^i    J       s+l  —  l'  '~s7.~ 

-  n  (r, , . .,r  )n, (r ) ]  6y(r)  dr  , 

which  since 

(11..9) 

may  be  written  in  the  more  convenient  form 


6£ 
s 


(Ii»-^i!s)  =  ^s+i^i'-'^sv^  -  na(£i'--'i,s)ni(£)] 

(I4..IO) 
[6Y(r)-(l/s)  ZI  6r(ri)J  dr. 
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To  take    advantage    of   the   asymptotic   expression   (2.22)    or 
(3.23),    (I4..IO)    is   further   decomposed  as 

-a,  ^  n%    2  ^s(Sl>-''£s)f>nl(r)  .. 

6ns(£i"-'£s)  +  ifn  3^ J-gs— * 

x>[6Y(r)    -.JZMr^Jdp     =    j[ns+1  (£!>•• ,£g,r)   - 

-  V'l'"'*.'   nl(£}   +  Ip  n     3n ?n— ]    * 

v[6r(r)    ~  7  ZI  ^(^i^J    dr    .  (1^.11) 

According  to    (3.1)    and   (3,12),    the  mean  compressibility 
)C     and  reciprocal   energy      (3     need  not  be    specified  more 
precisely,    for   the    combination 

nX/P     =     N(-n2  d2  In  Q/^n2)"1    ,  n  =  N/V   ,       (1|..12) 


A 


is   directly  determined  by      Q. 

By  virtue   of   the   relation    (3.23),    the    integral   on  the 
right   hand   side    of    (I4..II)    has   contributions   only  from     r 
within   some   conservative  multiple    of    the    correlation  length 
from   at   least   one    of     £,  ,..,r   ,      certainly  a  finite   region 
R     independent   of      V     when  all     r,      are   far   from  the  walls. 
Hence    (i|..ll)    is    seen  at   once   to  be    a  local  expression  for 
the    change   of     n    .      Further   simplification  may  be   achieved. 
Assume  that   the    spatial   rate   of   change    of      ^Y^f)     roay  be 

garded  as  uniform  in  region  R.   Then  6y(r)-(l/s)  ^~6y(rJ ) 
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re 


may  be  replaced  by   (r-(l/s)  5[~  r4  )•  S?6r(r ) .   where   r 
is  the  center  of  mass  of  r,,..,r  .   Similarly,  assume 
that   n^(£J   a^d  ^n-^rO/Jn  are  essentially  constant 
in  this  region,  for  we  are  concerned  with  large  scale 
inhomogeneities,  while   ns+1   is  also  uniform  in  that  it 
depends  only  upon  interparticle  distances.   Under  these 
conditions,  if  one  makes  the  center  of  mass  reflection 
r  — >    (2/s)  >  r.  -  r,   the  second  bracket  in  the  integral 
reverses  sign.   On  the  other  hand,  for   s  =  1  or  2,   a 
center  of  mass  reflection  of  r  leaves  the  distances  of 
n  ±ilrn,..,r  ,r)   unaltered,  so  that  the  first  bracket 
is  unchanged.   Hence  the  right  hand  side  of  (ii.ll) 
vanishes,  and  we  are  left  with 

(Jx.13) 
"  7  Z  6Y(r1)Jdr)  ^ns(r1,..,PB)/^n 

for   s  =  1  or  2.   This  is  of  0(1)   with  respect  to   N. 

Consider  first  the  case   s  =  1.   Eq.  (Ji.13)  then 

reduces  to 

V   o  ,i  ^n,  (r )  ^n,  (r,  ) 

5ni<£i>  "  §|  * 2J-3iT-"*<5i>  "  <*<£>>*  -jr^"  "MiW 

which,  noting  that   (<J/dn)  J  n^r )  dr  =  ^N/^n  =  V  and 
inserting  (2.20),  may  be  solved  for   6y(r1)   in  the  form 
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0        <^n_  (r,  )  -,  „  ^nn  (r ) 

*(&>     -    -|=  3^-feVjH  +IHF  <*<r)  dr 

Thus,    the   change   in     y(r,  )      required   to   produce   a  given 
change   in     n, (r, )      is,    except   for   a    constant,    identical 
with  that  which  would  be    obtained  by   changing      pu-     by 
means   of   a   change   in  total  particle   number    (or  overall 
density     n)      adjusted   to   yield  the    actual   local  density 
n,  (r, ) .      This   shows   that   a   slowly  varying   density   can  be 
interpreted  by   representing  each  fluid  element   as   an  open 
system  exchanging  particles  with  neighboring  fluid  elements 
(total  particle   number  being  maintained),    and  gives  meaning 
to      (l/6)y(r)      as   a  local   chemical   potential,    up  to   a 
constant   independent   of     r      (but  which  may  depend  on   the 
function     y) . 

Let   us   illustrate   this   by  considering   the   variation 
of  density  of   an  ideal  gas    in  a  uniform  gravitational 
field  in  the    z-direction   (a  potential   of     mgz )      due    to    a 
change    in  the    gravitational   force   constant     g.      According 
to    (I4.. 5),    this  will  result   from     6y(r)    =   -    p  m  6g  z. 
Further,    since      n(r)    =  Npmg  exp(-pmgz)      for   a  unit   area 
of  gas   between     z  =   0      and     z   =  00  ,      then 
6n(r»)    =    (l/g-pmz-,  )n(r,  )6g,      which  would  equivalently  be 
produced  by   a      6N      of      (l/g-pmz,)   n   6g. .      Now  the    chemical 
potential   is   given  here   by 
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[i  =  Jf/Jx  =    (<J/^N)(-l/p)[|  N  In    (2ron/ph2)    -  In  Nl    + 

exp  (-(3mgz)dz] 
=  const.  -  (1/0)  In  N  +  (l/p2mg)   , 


+  N  In 


r*> 


P       P 

so  that   6(Pn)  =  6N/N  -  6g/(3  mg  .   Under  the  above 
changes,  this  becomes   6(p|x)  =  -  (3m  6g  z.. ,   which  is 
indeed  identical  with  Sy^-i)* 

To  complete  the  picture,  consider  next  the  case 
s  =  2,   in  which  our  principal  interest  lies.   Observing 
that  ^T]  6y(r.  )  =  s  6  y(r)   through  first  order  gradient 
terms,  we  find  that  Eq.  (I4..I3)  takes  on  precisely  the 
same  form  as  Eq.  (I4..  llj.)  for   n-,(r).   It  follows  then  that 

.  c)rL(r,  ,r0)   <}n,  (r ) 

A 

In  other  words,  the  change  in  np(r-,,rp)   produced  by  6y 

is  equal  to  that  which  would  be  obtained  by  only  altering 
the  number   of  particles  in  the  system  so  that  the  local 
density  n,  (r)   at  the  center  of  mass  r  =  "x   (r,+r0) 
attains  its  actual  value.   Thus  we  reach  the  desired 
conclusion  that  a  variation  of  n, (r)   also  has  the  same 
local  effect  on  n?   as  an  equivalent  change  in  local 
density  achieved  by  changing  the  total  particle  number. 
The  infinitesimal  process  leading  to  (lj..lij.)  can  be 
iterated  as  long  as  neither  the  gradient  of  the  effective 
external  potential   y(r)   nor  the  form  of  the  two -body 
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distribution  vary  significantly  within  a  correlation 
length.   Starting  with  a  uniform  motionless  fluid  (of 
small  [surface  *  correlation  length/volume]  ratio),  we 
can  further  conclude  in  this  fashion  that  with  the  sole 
condition  that  the  external  force  and  local  velocity 
(and  hence  local  density)  be  slowly  varying  in  space, 
the  2-body  distribution  will  be  precisely  that  of  a 
uniform  fluid  of  density  equal  to  the  local  density  at 
the  center  of  the  two  particles. 

The  simple  expression  (I4.. II4. )  is  not  valid  for  nQ 
with  s  >  2.   Indeed,  the  vanishing  of  the  short-range 
contribution  to  the  right  hand  side  of  (I4..II)  was 
established  by  an  elaboration  of  the  familiar  argument 
[13]  that  for  n0(r,,r0)   there  exists  no  vector  of  local 
character,  symmetric  in  r, ,r_  which  can  combine  with 
Vy(r)   to  give  a  first  order  correction.   This  argument 
fails  for   s  >  2.   Instead,  a  density  gradient  can  provoke 
an  anisotropy  in  the  local  s-body  distribution.   Replacing 
(I4..13)  by  (I4..II)  in  (i;.  16),  and  representing  ns+]/ns  bY 
the  superposition  approximation  [II4J,  results  in  a 
relatively  simple  correction  to  (I|..l6)  which  we  cite  for 
the  sake  of  completeness:   if 

*s(r/L'">£s)~ri  n^r^  n^^n^r^r.Vn^r^n^r.)), 
then 
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dn   (r-.  , .  .,r    )     ^n,  (r  ) 

5n(r,,..,rJ    -    (     »  "V =S-/ Ip.)    5nn  (?) 

s  *w.X*        ^s  <^n  <^n  l  ~* 


n0(r, ,r ) 


(li.17) 


=  K    (£i»"«£b,J  [n(njr~)n~(r 


)    - 


n0  oo  (r,  ,r)  _  _ 

n,  Ir.  Jn,  lr;  1  ~       -,  ~       ~  -» 

A  A 

Here     n  is    the    asymptotic  form   of     n?. 

We   proceed  next   to    the  case   of   a   slowly   varying  local 
temperature,    anticipating   that   low  order  distributions 
will   again  maintain  their   equilibrium  form.      Changing  the 
function      p(rj      in   (li.lt)    has   two   effects   on     ng      of    (lt.6). 
First,    there   are    one-body  terms   which  will   alter,    namely 
y(r)      and  the      exp[-  4  ln(27imkT(r )  )  J      arising  from  normali- 
zation  of  the  momentum   distribution.      We  have    just    seen 

A 

that    a   slow   alteration   of    a   one-body  term   leads   to    an     n~ 
in  which  this    term  may  be    taken   as    literally   constant. 
There   remains   then  the   effect  upon  the   two-body  potential 
terms,    that    is,    upon 

p(q,r)   0(q,r)      ■     \   (p(q)+p(r))   jZKq,xO    . 

Clearly 

6(Bfif)      =      i    (6p(q)    +    6p(r))    0(q,r)    . 
d.  ^  -^  -»  /-» 

It    is    now   only   necessary   to   extend  Eqs.    (lt.8  )-(l|.10)    to 

two-body  variations.      Doing   so,    we   have 


-  38   - 


5*- S*±»  •  •  '£s  ^ 


/\ 


2   3  ~ij/ r^-  =  n  (r-.,..,r  )  ^1,  /,  6(r,-q)  6(r  -r)  (ij..l8) 

6pj3(q,r)        s  ~1»   '— s  *—  If 2       —1  ~        ~J  *" 

+  ^^l'-'^s'S'i^s^l'*-'^  *2(i'Z}  ' 
from  which 

26^^, .  .  ,ra)  =  I6p0(q,r)[ns+2(r1,..,rs,q,r)  - 

-  ^s(j:i»"»£s)S2(i'5:)j  d3  d£  +  2 j  ^  5P^(ri'a)[^s+itp-»rS'S)" 

-  SB(p1,..#rB)  ni(q)J  dq  +  Sn^r^..,^),  (^19) 

*  J  z  6P%i^'ni(a)ds +  ^s(ri^-^s)  SI'  °p^ti»sj)  » 

where  ^Z   denotes  the  omission  of  i  =  j. 

Our  inquiry  now  concerns  the  extent  to  which  the  change 
in  n  (r,,..,r  )      differs  from  that  which  would  occur  due  to 
a  uniform  change  in  temperature  corresponding  to  that  at  the 
midpoint  r  =  (l/s)  5~  r .   of  r,,..,r  .   Thus,  being  more 
explicit,  we  wish  to  consider 

nMr,,..,r  :p(r)+6p(r))  -  n  (r,, . . ,r  ;p(r )  +  6p)  , 

the  difference  of  two  types  of  variation,  where   6p   is 
the  uniform  variation  6p(r),   and  Y(r)   is  the  same  for 
both  distributions.   Again,  this  difference  may  be  divided 
into  an  asymptotic  and  a  residual  part: 
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(n3^/Np)^ns(£l,..,rs)/an[  J^(q,r,r)n2(q,r)dq  dr  + 
+  2f  ZA(q^±.7)  ^(q)  dq]  +  T]L  +  T2  +  T3  +  T^  , 
where 

+   (n3X/Np)  ^n    (r,  ,..,r    )/9n  <3n5(q,r  )/<)n]    dq  dr 

T2  •  J  ZA(q,P±,p)    CnB  (riffs')   Vi3    "  *s+l(rl'*-'iVS;)  + 
+    (n3^/Np)   t)n   (r,,..,r    )/^n  Jn,  (q)/c)n]    dq 

T3  =     "  ^^'"^s'l^^S'^  nl(i?   di 

T,    =     -  n(r,,..,r    )   Jl' A  (r,  ,r  .,?)  (lj..20) 

and 

A(q,r,?)      =      [i   fip(q)    +^   6p(r)    -    6p(r)J   0(q,r)    . 

The    "unperturbed"   distributions   all  have   temperature 
dependence      p(r),      although  this   is  not   explicitly 
indicated.      The   elimination  of   correction  terms     ti"tk 
in    (I4..20)    is   possible   once  more   only   if      s   =  1   or   2, 
which  we  henceforth  assume. 

In  order    to   reduce    (l|..20),    let   us   more    carefully 
define    the    slow  variation  of      6(3,       in  the    sense    that   if 
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L  is  a  macroscopic  length,  Jt     a  microscopic  one,  then 
LV6p   is  of  order   6(3,   but   XV6p   and   £2<7V6p  may  be 
taken  as  zero.   Consider  first  the  term  T. .   Again  assuming 
local  dependence  of  n q+p,  n?,  0     on  interparticle  distance 
alone,  both  p(q,v)      and   [n  rip-n  +~...J   are  unchanged  by 
an  inversion  through  the  center  of  mass:   x  — »•   2r-x. 
But  $(q,r)   insists  that   q  be  close  to   r   and  the 
asymptotic  deviation   [ri  n9-n  .  0» .  .  J   then  requires  that   q 
and  r  be  within  the  order  of  a  correlation  length  from 
|r,,r2y   Since  the  coefficient  of  $(q,r),  averaged  with 
its  inversion,  is 

|  [|  6p(q)  +  |  6p(r)  -  6p(r)J  +  |  [j   6p(2r-q)  + 

+  |  6p(2r"-r)  -  p(r))]  , 

with  leading  term  now  going  as   X  77p(r),   the  term  T, 
may  be  dropped.   The  elimination  of   T~   and  T-,   proceeds 
in  precisely  the  same  fashion,  while  for   T,   an  average 
over  inversions  is  both  impossible  and  unnecessary. 

Thus  only  the  asymptotic  contributions  to  (I4..2O)  remain, 
and  the  second  of  these  vanishes  for  the  same  reason  as  that 
ascribed  to   T-,.   Since  the  resulting  relation  is  equally 
applicable  to  variations  of  n0(r, ,r~)   and  n,  (r),   it 
follows  that,  as  in  (1^.16), 
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n2{Tlfr2;p{r)+&p(r))    -  n^r-^r^ptr  )+6p) 

dnp(rn ,rpJ  dn, (r) 
=   (    gn    /   in^  )^1(r;p(r)4-6P(r))-n1(r;p(r)+6p)J< 

Combining  with  ([(..16),  we  see  that  a  uniform  change  of 
reciprocal  temperature   6p  followed  by  an  alteration  of 
y(r)   sufficient  to  bring  n^rj (3+6(3 )   to  its  true  value 
n^rjp+ep)   results  in  a  value  of  n2^l»r2^   identical 
with  that  produced  by  the  non-uniform  change   6(3 (r).   This 
infinitesimal  process  can  then  be  iterated  to  show  that  a 
slowly  varying  temperature,  with  no  restriction  on  the 
magnitude  of  its  change  in  the  large,  can  be  reduced 
insofar  as   np(r,,r2)   is  concerned  to  a  uniform  terperature 
p(r)   with  modified  y{r)      and  thus  to  the  corresponding 
equilibrium  distribution  at   p(r),   n,(r).   As  a  consequence 
of  the  analyses  of  (J+.16)  and  (I4..2D,  the  validity  of  (I4..7) 
for  s  =  2  has  now  been  demonstrated. 

In  Appendix  D,  we  show  that  the  pattern  of  inference 
established  above  may  be  reversed,  in  that  the  asymptotic 
forms  (2.22,  3.23)  are  themselves  consequences  of  the  local 
character  of  distribution  functions  together  with  the  basic 
Ornstein-Zernike  relation. 
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5.   Conclusion. 

We  have  shown  in  this  paper  how  the  assumption  of  the 
existence  of  a  finite  correlation  length  in  a  fluid  yields 
explicit  expressions  for  the   (l/N)   terms  in  the  joint 
distribution  of  two  sets  of  particles  which  are  far  apart 
compared  to  the  correlation  length.   The  form  of  these 
terms  was  then  utilized  to  prove  the  local  nature  of  the 
low  order  distributions  in  a  system  with  spatially  varying 
intensive  parameters. 

It  appears  to  us  that  the  central  problem  in  the  theory 
of  equilibrium  fluids  is  the  proof  of  the  existence  of  such 
a  length.   As  mentioned  in  the  introduction,  this  is  related 
to  the  distinction  between  fluids  and  crystals  and  hence  to 
phase  transitions  between  these  forms. 
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Appendix  A 

In  this  appendix  we  attempt  to  give  a  more  rigorous 
mathematical  formulation  of  the  concepts  of  finite  correla- 
tion length  and  asymptotic  form  of  distribution  functions. 
It  is  clear  that  these  concepts  can  be  made  precise  only 
when  we  have  some  procedure  for  letting  the  particle 
number  N  approach  infinity.   When  the  system  is  completely 
uniform,  i.e.  periodic  boundary  conditions,  this  passage  to 
the  limit  is  indeed  straightforward:   N  -»■  oo  ,   V  -*>  co , 
N/V  =  n   (although  even  here  one  may  have  to  specify  the 
ratios  of  the  various  sides). 

For  a  non-uniform  system,  we  may  imagine  the  restric- 
tion to  successively  larger  volumes   V  by  imposition  of  an 
appropriate  (short  range  and  infinite)  wall  potential.   All 
other  conditions  are  to  be  held  fixed  at  predetermined  values 
throughout  space:   internal  and  external  (e.g.  gravitational) 
potentials,  as  well  as  local  (for  a  local-equilibrium  ensemble) 
temperature  and  chemical  potential.   We  assume  then  that  these 
quantities  are  sufficiently  bounded  that  if  the  particle 
number  N(V)   for  each  volume  is  suitably  chosen,  the 
distributions  n.  (r, ,  .  .,r,  ;N)   of  finite  order  k  will 
possess  limiting  values  n,{r-.f.,,r^)      as   N  -*-  oo  .   The 
value  of  nk(r1}..,r,  )   may  of  course  depend  on  the  way 
the  volume  becomes  infinite.   In  the  following  we  shall 
think  particularly  of  two  ways  of  passing  to  the  limit: 
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the   first   is    such  that   the   positions   of   all  k-particles 
become    infinitely  far  from  the  walls,    the   second  that    some 
of  the  walls   remain  located  at   a  finite   distance   from  the 
k-particles.      As   a  prototype   of   the   first    situation  we    start 
with  a   cube   defined  by      -  h  L  <  x,y,z   <  ^  L     and   then  let 
L   ->   oo ,      while   the   second  may  be   represented  by   starting 

instead  with  a   cube      0  <  x  <  L,      -2L-y»z-2L*      The 
second   situation  is  necessary   for   the   examination  of  wall 
effects   on  the    distributions   in  the    limit   of   infinite      L 
for   an  otherwise   uniform  system. 

We  will  now   specialize   our   discussion  to   the    one   and 
two-particle   distributions.      The    two-particle   Ursell 
function  for   given     N,V     is 

4 


'0(r,,r0;N)   =  n~(r,,r~;N)  -  n(r-.  ;N)n(rP;N) .  (A.l) 

d      — 1' — d*  d    S^L'^d'  ~.±'        ~d 

P  vanishes  when  either  of  the  particle  coordinates  is 
outside   V  and  has  the  property  that 

Jv  O^    (r1,r2;N)  dr2  =  -  nCr-^N)  .     (A. 2) 

Our  previous  assumption  assures  the  existence  of 

N^oo  2~l'-2'       ^2~1'~2  (A#3) 

=  n2(^l'i:2)  "  n(^l}  n(£-2}  * 
The  Ursell  function  v/.(r,,r?)   is  defined  in  such  a 
fashion  that  it  vanishes  for  statistically  independent 
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particles.   The  correlation  length  of  a  fluid  is  therefore 
related  to  the  scale  on  which  */"„   vanishes  as   lrT?l 
increases.   The  possibility  of  correlation  over  large 
distances  may  be  intrinsic,  as  in  a  crystal  or  at  the 
critical  point,  or  may  be  due  simply  to  the  constraint 
of  a  fixed  number  of  particles  in  the  system.   It  is 
clearly  the  second  part  of  this  long  range  correlation 
with  which  we  are  concerned  in  this  paper.   The  magnitude 
of  this  non-intrinsic  correlation  will  vanish  as  the  size 
of  the  system  increases.   Hence  in  the  absence  of  intrinsic 
long  range  correlation,  which  we  here  assume,  we  will  in 
the  limit  have  for  some   m 

J  r^|J^2(r1,r2)  |  dr2  =  0(1)   for  0  <  A  <   m.  (A.lj.) 

The    value    of      m      determines    the   rate    at  which    J  p 

approaches    zero   as      r-,2   -*■  oo  ;      if    (A.ij.)    holds   for   all    X, 

the    approach  is    exponentially  fast. 

In  addition   to    vanishing  rapidly   at   large    separation, 

>-/  »(r,,r0)      will   onlv  have    a   small   fractional   deviation 
2  ~l'"-2 

from     ^/"0(r,  ,r0;N)      at   small     r-,0.      Hence   the    difference 

^2(£l^2^N)      ~    ^2(£l^2'N)    "    ^2(£l>£2}  (A'5) 

(defined   as    vanishing   outside   of      V)      may  be    termed   the 


asymptotic   part   of     >^"_(N).      In   applications,    use   is   made 


of    the    short-range    character   of     <Zjp,      coupled  with 
properties    of    the    integrals   of     0?(N)      over    the    full   volume, 
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or  more  precisely  of 

lira   /  0L(r,  ,rQ;N)  X  (**.,)  dro   >  (A«6) 

where   /(r?)   is  a  bounded  function:   |/£(i*p)|  <  M. 

We  shall  first  discuss  the  value  of  the  integral 
(A. 6)  when  /C(r0)      is  a  constant.   This  will  yield 
sufficient  information  about   0~(N)   that  the  properties 
of  (A. 6)  will  be  obtained  with  few  further  assumptions. 
The  advantage  of  the  restriction  to  X-      being  a  constant 
is  that  since 

h*Z    (£l'£2jN)  ^2   =   ■  n(£l*N)  "  Jv^(£l^2)d£2   (A*7) 

(employing    (A. 2)    and    (A.5)),    then   only   the   integral    of     <•„ 
is   required,    which  being   short-range    can  be   found  by 
integrating   over   a   finite   volume   as      V  -*■   co  . 

Consider  now  the   process   discussed  in  Section  2   of 
dividing   the   volume      V     into     VA      and     VB  =   v  "   VA- 
Using  the   previous  notation,   we   have 

^2(£l'£2'N)      =    ^n2(£l>£2lNA'N~NA> 

(A. 8) 

"  <n(£l>  lNA»N-NA><n(i'2)  'NA>N-NA>  ' 

If     ^7^(r,  j^TpjN)      is    integrated   over      VB     with     r-,      in 

VA,      then  from    (A. 8) 
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/vB^l'^N)d£2      ■  <'(N-NA)n(r1!NA,N-NA)> 


"  ^-NaX^i^a'^a^  (a-9) 

=   -  <NA  n(r1|NA,N-NA)>  +  <NA><h(£l  |NA,N-NA)  > 
oo    <(6N    )k>    nk  . 

k=2      (NA)ic  X     KI        ^n*     A        X     A  A 

where      nA   =  NA/VA    ,       6NA  =  NA    -   NA    . 

By  taking  the    limit   of      N,V  -*>   oo  ,      with      VA     remaining  fixed, 
it   follows   that 

f  <T0    (r-,,10    drP     =      f       11m   J^p(r,,rpjN)    drp  (A. 10) 

VA  A 

=   lim        S^2^lfv2^)    dr2      =    -ntr^-lim   /  ^(r^jN )  dr2 

N-*oo       VA  N-»oo     Vg 

oo    <(&NA)k>nk        jk     _  _ 

=   -(£,)   +  ZZ  — fcj-  ^  [-^   [nAn(r1|NA,N-NA)]    - 
k=2    (NA;  ^nA 

_     ^kn(r]L|NA,N-NA) 

"   nA  73 J    • 

dnA 

where  all  quantities  on  the  right  hand  side  are  to  be  taken 
in  the  limit  of   V  becoming  infinite. 

Next  let  N.   increase  to  infinity.   The  quantities 
inside  the  square  brackets  in  (A. 10)  are  of  at  most   0(1), 
and  this  remains  true  if  for  the  purposes  of  rigor  the  series 
is  truncated,  the  final  term  being  evaluated  at  other  than  NA . 
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We  now  assume  that  the  fluctuations  in  NA   are  bounded 
in  order  of  magnitude  by  the  fluctuations  in  an  ideal  gas, 
where  it  can  be  shown  that 

<(6NA)k>  <   0(NAk/2)  .  (A. 11) 

This  yields,  for   V.   approaching  infinity  in  such  a  way 
that  r,   becomes  infinitely  far  from  the  interface   SA 
(although  not  necessarily  far  from  the  container  walls), 

f^(r,,rJ  dr„  =   -n(r,  )+n2kT^  ^n(r,  )/<)n  +  o(l )     (A. 12) 
VA 
where  we  have  defined 

<(5N.)2> 
lim     — —^ =     n  k  T/  (A. 13) 

NA-*>oo  NA 

with  the    understanding  that   to   order      o(l)      in     NA, 
nkTX(NA)      and     n(r1,NA)      coincide   with  their   limiting 
values.      But   according  to    (A.I4J, 


f^l^2]    dZ2      =   f^Z^l'^    dX2   +   0(1)    '  {A*1^) 

A 

It   follows   from    (A. 7)    and    (A. 12)   that 

lim     /  ^2(r1,r2;N)    dr2     =      -  n2kT/£  ^n(r1)/^n   ,  (A. 15) 

N-*»oo        ""*      **  "* 

which  is  the  result  desired. 

At  fixed  r,  ,   characterization  of  #L(r, ,r0;N) 
will  entail  two  regions  of  interest.   For  v~      separated 
microscopically  from  r, ,   the  general  argument  following 

-  I4.9  - 


(2.10)  suggests  that  the  effect  of  the  boundary  may  cause 

,5^(N)   to  deviate  from  J^  by  the  order  of   (l/N)0(N2/^)  = 

-1/3 
=  0(N  '~>).      For   r2   outside  an  intrinsic  correlation 

length  from  r^,   the  effect  of  the  constraint  (of  precisely 

N  particles  in   V)   could  still  introduce  terms  of  order 

-1/3 
0(N  '  ).   However,  this  could  certainly  not  be  true  in  a 

uniform  sense,  due  to  the  existence  of  the  integral  (A.15). 

We  would  like  to  argue  that  in  this  region,  the  dominant 

term  in  j?L(N)   ^s  given  by  the  second  term  on  the  right 

hand  side  of  (2.19),  which  is  of  order  0(l/N).   This  is 

certainly  compatible  with  Eq.  (A.15)  and  as  we  shall  see 

is  almost  implied  by  it. 

A  very  strong  argument  for  this  behavior  of  JZL(N) 

can  be  made  by  considering  rather  than  j2L(r,  ,r~;N)   a  very 

similar  quantity 

^l'^V   E  ^2(^1^2|IA)  -^L2(^1^2) 

(A. 16) 

where     3~Av~ ,  r„  |N.  )      =      lim    3~0(r,  ,r5  |N  A,N-N.  )    . 
2  ~±  ~d      A  n-^oo       d   '  L    'd      A 

j^p(r,  ,r„  |N.  )      is    the    deviation  of  the    Ursell   function  for 
an   infinite    system  without    any   constraints    from  the    one    for 
the    same    system  when   a  given   volume      V.       is    constrained   to 
contain     N.      particles.      As   far   as   the    constraint   of   a  fixed 
number   of   particles    is    concerned,    the    subsystem   inside      V. 
may  be    regarded   as    a   closed    system  with   a   special   type   of 
boundary.      The   function    J'JrprJI^      should  therefore 
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mimic  for  r, ,r~   inside   V.   the  behavior  of  the  function 
jTpCr-,  ,rp;N.  )   for  a  closed  system  of  N.   particles  inside 
the  container   V..   Indeed,  as  far  as  the  integral  over   V"A 
is  concerned,  the  two  functions  have  the  same  behavior: 

A 
The  properties  of  0p^£l*£2^A^   can  be  found  from  an 
expansion  similar  to  that  leading  to  (A.9),  which  yields, 
as  a  development  in   (l/N. ), 

4^2'V   ■  A(-Sl^2)-<(<niA)^t^*?2(Sl*^A) 

^^A 

+  -^-  n(rj  |N,  )  -^-  n(rJI,)]  +  ...  . 
*\  ^A   " 

For  rn   and  r„  far  from  each  other,  this  expression 

reduces  to 


(A. 18) 


NA        *nA    0nA 


if  r-.   and  rp   are  both  interior  to   V.   (the  correction 
term  vanishes  otherwise).   The  right  hand  side  of  the 
expression  is  indeed  of  the  same  form  as  that  found  in  (2.19) 
for  the  function  #L(r-,  ,r2;N.  ) . 

Let  us  now  return  to  the  discussion  of  #L(r..  ,rp;N) . 
On  the  basis  of  (A. 15),  we  divide  0p(N)   into  two  parts: 

02(£l»£2'N)   =  N_1^20)(^1^2;N)  +  ^21)(-1'£2'N)  '       (A. 20) 
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Here      #L      (N)      is   to   remain  bounded  as      N   -*-   oo  ,      i.e. 
^0)(N)    =   0(1),    and   satisfies 


i     0  Acs\  p      -Jn(r,  ) 

lim       ±  J   ^2      (£l'52'N)    d^2     =      "n  kT/t!      ^      ;  (A. 21) 


the   remainder     j2L      (N)      may  be    of  order  greater  than 
0(l/N)      but  by   virtue   of    (A. 15)   must   satisfy 

lim      j   0^1Ur1,ro;N)    dr0     =     0    .  (A. 22) 

N-*co         2       ^  ~2  ~2 

The  decomposition  (A. 20)  is  of  course  not  unique,  since 

one  can  for  example  add  and  subtract  any  term  of  order 

o ( 1/N ) . 

The  vanishing  in  the  limit  of  the  integral  (A. 22) 

can  occur  in  two  basic  ways.   First,  the  part  of  0„      (N) 

which  is  of  order  greater  than   0(l/N)   may  be  effectively 

restricted  to  a  finite  region.   Second,   JZL   (N)   may  be 

of  order  greater  than  0(l/N)   over  the  full  volume,  but 

have  an  oscillatory  character  leading  to  a  vanishing 

integral.   We  assume  now  that  in  a  fluid  the  latter 

situation  does  not  arise  and  express  this  formally  by  the 

condition  that  a  decomposition  (A. 20)  exists  for  which 

lim  f\0^hrlfv2^)\    dr2  =   0  .  (A. 23) 

N-»oo 

This  assumption  implies  that 


lim  )X(rltr2)   ^(r-^r^N)  dr2 
N-*-oo 

r(0) 


(A.214.) 


«»   lfr(?1^2)  *2   (£l'~r2'N)  *S2 

N-*oo 
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where      /-      is   an   arbitrary  bounded  function  of    _r,      and     r~. 
Hence    in  all   applications   in  this   paper  we   need  only  be 
concerned  with  the   value    of     ftp      (r-,,r2;N). 

Consider   now   a  homogeneous    system.      Let  us   denote   the 
limiting  value    of     02O)(N)      by     $^    , 

lira       j40)(r,,rPjN)      =     ^0)(r-.,r9)    ,  (A. 25) 

N-»oo      2       ~L  ~d  d       ~X  ^d 

and   its    asymptotic  part  by     \|/    , 

11m     ^0)(r,,r~)      =     t(r,  )    ,  (A. 26) 

(p2I-oo    2       -1-2  -1 

both  of  which  limits  we  assume  to  exist.   The  existence 

of  these  limits  means  that  we  discard  the  possibility  of 

£L   (N)   having  an  oscillatory  behavior  extending  the 

length  of  the  container.   The  limit  \r~  \    -*■   co   is  here  to 

be  taken  in  such  a  way  that   r~  moves  infinitely  far  from 

the  walls.   Combining  Eqs.  (A. 21)  and  (A. 26)  yields 

^(r^      =   -n3kT/  ^nd^J/^n  ,  (A. 27) 

so  that   finally,    using    (A.2l|), 

I^(Sl»52)^2(^l'£2'N)    ^2     =/^(£l'52)^2(£l'J2)   d^2        (A-28) 

-  n3kT  X(^n(r1)/^n)[N"1(f^(r1,r2)dr2J    +   o(l)    . 

For  a  system  which  is  not  completely  uniform,  the  limit 
(A. 26)  may  still  exist  if  the  system  is  asymptotically 
uniform,  in  which  case  (A. 28)  clearly  remains  valid.  However, 
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if  the  region   | r~ 1  -»>  co  is  not  uniform,  (A. 26)  will  not 
hold  and  must  be  replaced  by  the  strong  assumption  of 
statistical  independence  of  distant  parts,  i.e. 

lim   j40)(r,,rJ   =  ♦(r,)  ¥(xO  .       (A. 26') 
|r12 1  ->oo  2   -1  ~2       -1    ~2 

This  is  perhaps  the  strongest  assumption  we  have  made  and 
it  may  possibly  be  violated  at  low  temperatures.   If  we  do 
make  this  assumption,  then  (A. 28)  again  holds,  and  we 

further  have  on  combining  with  previous  equations 

,       -      _ ^n,  (rn  )  Jn,  (r0)  , 

fzW     -     §  AtfJ^i -2^*-  *  02   <N)   ♦  o(l/N)  (a>29) 

where       f\0'2(N)\    dr2  =  o(l). 

The  generalization  of  the  above  discussion  to  higher 
Ursell  functions  is  quite  direct  and  will  not  be  carried 
out,  but  the  final  expression  will  be  written  down  in 
Appendix  C. 
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APPENDIX  B 

In  this  appendix,  we  extend  Eq.  (2.22)  for  the 

asymptotic  value  of  a  distribution  function  to  fluid 

mixtures.   This  result  is  then  applied  along  the  lines 

of  Ornstein  and  Zernike  to  the  scattering  of  visible  light 

from  mixtures.  For  simplicity,  we  shall  consider  explicitly 

only  the  case  of  a  two-component  fluid  consisting  of  N 

atoms  of  one  kind  and  N,   of  another  kind. 

b 


The  joint  distribution  of  k,   particles  of  species   a 

2 


and  k0   of  species  b   will  be  written  as   n,    .   The 
2  KjK. 


asymptotic  value  of  n  +ff         +0       when  the  set  consisting 

ql  Zl   q2  XZ 
of  q-,   particles  of   a  and  q~  particles  of   b   is  very 

far  from  the  other  particles  can  be  found  by  an  investigation 

similar   to  that  used  for   the    one-component   fluid.      It   yields 

VA  VA  *  V2  n4V  If  <kT  ^'^V  ~?V 

where      a     and      p      can  assume    the    values      a      and     b,      and 

^a/fop     =     ^Np/^^a   '  (B*2) 

When  a  parallel   incident  beam   of  radiation  of  intensity 
I        and  frequency     v     falls   upon   this   fluid,    the   intensity 
of   the   radiation  at   distance     r     which  has  been  scattered  at 
angle      9     may  be    computed   in  a  way  entirely   analogous   to    that 
used   for   a   single    component   fluid. [15] •      It   is   given  by 
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1(9)      =     i   (1  +  cos2  9)   I      (r  /r)2  x 

d.  0  0 

(B.3) 

{Nafa+Vb  +W[ranaa+2Vbnab+fb'W   ^"^  *i  <*2   • 

Here  f   and   f.   are  the  atomic  scattering  factors  for 

single  atoms  of  type   a   and  b  respectively,   r  =  e  /m  c2 

o       e 

is  the  classical  radius  of  the  electron,  the  n  „(r, .r«) 

a(3  — 1'~2 

denote  the  two  self-  and  one  mutual-  two  body  distributions, 
and  K  is  the  change  in  wave  vector  on  scattering 
(K  =  (l^v/c)  sin  |  0). 

In  the  above,  we  can  subtract  a  constant  term  from  any 

na«   since  for  a  system  of  volume   V  this  will  contribute 

-1/3 

only  when  K  ^V  '  -1      and  is  thus  indistinguishable  from 

the  transmitted  beam  as   V  ->  oo  .   We  shall  therefore 
subtract  the  corresponding  asymptotic  value.   The  integrand 
in  (A. 3)  then  exists  only  for  small  values  of  r,p,   of  the 
order  of  the  correlation  length.   Further,  if  the  system  is 
uniform,  then 

nap(^1^2}   =  na  np  Sap(r12)  ,  (B.l*) 

and  Eq.    (B.3)    can  be  rewritten  as 

I(s)      =     |   (1  +   cos2  0)    Io(ro/r)2  V   j(s)  (B.5) 

where 

oo  >- 

J(s>   ■  Va  +  Vb  +    Lt^a^aa'-'-Saa100^ 

+  2'aVanbt8ab(r)-8ab('°>J+'fnbt8bb(«-)-8bb(»»}^2<5^»ar 

and      s   =    (Ipiv/c)    sin  -^  0  (B.6) 

-  56   - 


(reducing  correctly  to  the  one -component  case  when 

fa  =  f>      gap  =  g»   na+nb  =  n)*   Eqs*  (B'^J»  (B*6)  are 
appropriate  for  determining  the  long  wavelength  limit. 

If   v   is  sufficiently  low  that 

c/v  » JL  (B.7) 

where    £,    is  the    correlation  length   --  beyond  which  the 
integrand  in   (B.6)  must   vanish  --   then      s  «   l//C,      and  we 
may  replace      s   n  rs     by     1.      Writing 

nap(oo)      *     nanp   "    (1/V)    ^ap   >  (B'8) 

it  readily  follows   that 

j(s)    -*    j(0)    =   Yaa   t\  *  Ybb    4   *   2Yab    fa  fb    .  (B.9) 

But    according  to    (B.l),   we    can  write 

YaR      =     kT    (dn Jdv- s)|  •  (B.10) 

ttP  a        P    lT,V,txa 

Eqs.  (B.9)  and  (B. 10)  constitute  thr  result  desired. 
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APPENDIX  C 
We  now  extend  Eq.  (2.22,  3-23)    to  the  asymptotic  value 
of  the  intrinsic  correlation  functions.   These  may  be  defined 
by  the  sequence  of  relations 

J-^r)      =     n^r) 
v7"„(r,r«)  =  n5(r,r')  -  n,  (r)  nn(r') 

(c.D 

^T,(r,r',rn)  =  n,(r,r',r")  -  n3(r,r')  n,(r")  -  n0(r,r" )n. (r ' ) 
-  n9(r',r")  nn(r)  +  2nn  (r  )nn  (r  •  )n_  (r"  ) 

•   •   • 

The  fundamental  property  possessed  by  ^J^    (r-i>««>r  )   is  that 

of  vanishing  in  any  region  in  which  the  set   r, , . . ,r 

decomposes  into  two  or  more  independent  subsets,  i.e.  such 

that  the  distribution  functions   n.   correspondingly 

decompose  into  products.   This  is  because  all  correlations 

due  to  subsets  have  been  subtracted  in  forming  <«*/*   from 

s 

n  ,   leaving  only  the  "intrinsic"  s-body  correlation. 

The  connection  between  the   n   and  the  ^   may  be 

s  v   s  * 

defined  more    concisely.      For    this   purpose,    introduce   a 
"test    function"      f(r),      and   integrate    to   form   the    constants 
(functionals   of      f) 


n 
s 


[f]      =      /...I   n    (r, , ..,r    )    f(r.)..f(r)    dr,..dr 

<J  J        S^l'         *~S  —1  ~S  —1  ~£ 


(C.2) 
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and   the    combinations 


CO 


o 
oo 


nosl 


(0.3) 


The  distribution  functions  may  be  recovered  by  variational 
differentiation: 

-  rs 


as(£l— '£•>  "  5  n[fj/6f(r1)..6f(rs) 
^i**!'"'^  =  6S^[f]/6f(r1)..6f(re) 


f=o 
f=o 


(C.lj.) 


It  is  then  found  that  the  sequence  (C.l)  achieves  the  very 
concise  form 

J^[f  ]   =  In  n[f  ]  .  (C.5) 

The   required  extension  of    (2.22,    3.23)    is   now  readily 
accomplished.      From    (3.17),    we  have 


r^      _        (o)r--i        1  n%.  _2     J      _(o)r.Pi 
n[fj      =     n        If  J    -  tt  "5f  n     -*-*  n        LfJ    ...    , 


N   2(3 


(C.6) 


dn< 


and    (C.5)   becomes 


2  _(o) 


^/r^i      -     -i         (o)rfl               1          1  nX  _2  S   n"u'[? ) 
J[f]      =     In  n        [f]    -      (q)  y  jg  n  2 ... 


nvu/[fj 


<K 


=      ( 


l-^n2A)ln  n(o)[fJ 


N   2p 


(C.7) 


<K 


(o) 


1  n£  ,      d    In  nvu/[f]^2 


To   find      C?q+^      we  must   apply      5 q+~/6f  (£.,_).  .6f(rq+^  )  |f=0. 
But    suppose   that   the    sets      q,Jl     are   asymptotically   separated. 
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Then  according  to  (3.22),  the  n     decompose  accordingly, 
so  that  zh+j     obtained  from  In  n   [f]  must  vanish.   In 
the  same  fashion,  there  can  be  no  contribution  from  the 
[n  t)    In  n   [f]/^nj    term  unless  the  complete  set  of 
derivatives  corresponding  to  the  set   q  operates  on  one 
factor,  those  for  set  J-     on  the  other.   We  conclude  at 
once  that 

the  desired  result.   It  should  be  pointed  out  that  our  final 
expression  (C.8)  depends  only  upon  the  asymptotic  form  of 
the  distribution  functions,  and  its  validity  is  therefore 
not  restricted  to  the  gas  region  for  which  the  separation 
into  n     and  a  remainder  (necessary  for  (C.6))  is  possible. 
We  may  note  that  3"      becomes  a  higher  order  infinites!- 
mal  than  l/N   if  its  particles  divide  into  three  or  more 
groups.   Another  consequence  of  the  asymptotic  form  (C.8), 
which  was  previously  discussed  (see  Appendix  A)  for  the 
case  of  <Tp,      is  an  expression  for  the  integrals  of  ~q+. 
in  an  infinite  system.   Indicating  the  explicit  dependence 
of  Ursell  functions  on  particle  number  N,   we  have  generally 
(where   q  >  0) 

/4t«(5L-*«rt*',,)  "Vi  ••  *«*  (09) 

Consider   now   the    corresponding   integral    for   the    difference    of 
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0~  +     and  its  asymptotic  value   <J~a+^"    given  by  (C.8), 

^1  (as.  ) 

The  integrand  of  (C.10)  is  of  a  local  nature.   We  can 
therefore  pass  to  the  limit  of  N  -*•  co   inside  the 
integral  sign,  yielding  finally 


(CIO) 


/^  <£!.... £q+//»)    %+l   ••    <*q+* 


of  which  the    Ornstein-Zernike  relation  is   a   special   case. 


(C.ll) 
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APPENDIX  D 

Here,    we    "close   the   circle"    and   show,    using   the 
well-known  special   case      q  =   JL    =1,      that   the   asymptotic 
formula    (2.22,    3.23)    is    implied  by  the   local   character   of 
distribution  functions.      This  may  be   done   in  two   stages. 

First    consider   the    conditional   distribution 

n    (r,  ,..,r    |r    .-,  )      when     r    .,       is    far   from   the    set      r,  ,  .  .  ,r    , 
q  "1*       '~q'~q+l  ~q+l  <v>l*       '<*  q 

of   center   of  mass     F.      Now   the   effect   of   fixing     r  +,      is 
to   change    the    local    density  at     r      from     n(r)      to 


n«(r)      =     n(r|rq+1) 

=  n(?)    _  I        nX  n  frig)    jj^gg+lj 

Further,    if     n   (...:n)      denotes   the    -joint    distribution 
'  q  '  u 

function  for    a  uniform   density     n,      the    local  dependence 
of     n        on     r      tells   us    (see    (i|,  13,    [j..  17))    that 

for    some    appropriate   vector  function     A    .      Also, 

n   (rn...,r    |r    . -,  )      is   similarly  expressed,   with     n'(r) 
q  «-i*       '~.q'~q+±  ~*- 

replacing     n(r).      Hence 


(D.l) 


(D.2) 
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n   (r,  . .  .r    |r    .-  )    -  n   (r, ..r    ) 
q  ~1        ~xq'~q+l  q  ~±     *»q 

=  n    (r,  ..r    :n'  (r) )    -  n   (r-, .  .r    :n(r) ) 
q  *-.±     ~q7       ^^  q  *-±     a- q*     — 

+  Vn'(r)'A    (r, ..r    :n«(r))    -  *\7n(r)  •  A    (r,  .  .r    :n(r) )  +  . . 
v         ~        ^q  «v  1     *-q  "  »~        ~  q  *-^x     «*  q        -^ 

■    tn.(r)-n(r)]        1  ^^        =       /  ^  (D.3) 

+   [Vn'(F)-Vn(?)]'A    (r,..r   :n(r)) 

AA.   (r,  •  .r   ;n(r) )        iw/r:\ 
♦   [»•(£)    -  n(f»  Vn(r)  •    -<  S^       =       /  fcjgU.    . 

Since   the    only     r      dependence   of     n* (?)    -   n(l*)      is 

through  the   factor     ^n(r)/^n,      it   follows   that 

e)n(F)        ^n(F) 
Vtn'(r)    -   n(r)J   =    [n'(r)-n(r)J  V<~ g~"  /  -jjf-  +-    •      ^-^ 

Inserting   in   (D.3)    and  comparing  with     e)n    (r1..r    )/^n 
from    (D.2),    we   obtain 


V^l"1^^    "    nq(-rr^q) 


(D.5) 


d}n    (r,  ..r    !r    .,  )        d)n(r) 

Substituting  from     n' (r )    -  n(r)      of    (D.l)    and  multiplying 

by     n(r    . ■,  )      now  results   in  the      q,l      case    for   asymptotic 
~q+-L 

separation: 

Vlfe'^q-Vl1      ~*      VSl"*^   n(^q+l) 

1   nXn  ^nq(^l'^q}   n^n(i'q-H) 
N     (3  fljn  dn 

Next   consider     n»(r    .-...r    .  Jr,  .  .r    ) .      The   effect   of 

JL  ~q+l     »--q+/,~l     ^q 
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fixing  r,,..,r   is  to  change  the  local  density  at  r^ 
(center  of  mass  of  the  set  of  A-       particles)  from  n(r '  ) 
to   n  (r'lr, ..r  ),   which  we  have  just  obtained.   Applying 
the  above  argument  to  the  set  of  Jc      particles  then  recovers 
the  full  relation  (2.22,  3.23). 
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